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The relations between neutrino NC and CC structure functions and cross sections, which are driven by isospin 
symmetry, are discussed for nucleons and nuclei. 



1. Introduction 

The studies of neutrino interactions with nuclei 
has been an actively developing field both exper- 
imentally and theoreticaly |1I2| . Nuclear physics 
studies with (anti)neutrino beams are relevant in 
two different contexts. First, it must be noted 
that neutrino collision experiments usually use 
heavy-nucleus targets. For this reason a good un- 
derstanding of "standard" nuclear effects is nec- 
essary for interpretation of high precision experi- 
ments with neutrino beams. On the other hand, 
neutrino experiments can provide better under- 
standing of physics of nuclei that is important 
and interesting per se 

The interaction of (anti)neutrino with matter is 
mediated by charged or boson (charged 
current, CC), or by neutral Z boson (neutral cur- 
rent, NC). For this reason the studies of neutrino 
interactions provide information which is not ac- 
ccssable with charged- Icpton probes. For exam- 
ple, the neutrino NC interaction strength is de- 
termined by the weak mixing angle 9iv and the 
studies of relative rates of NC and CC neutrino 
reactions provide a tool for the measurement of 

In this paper we analyse relations between CC 
and NC neutrino structure functions for nucleons 
and nuclei. In Sect.[21we discuss the derivation of 
NC-CC relations for structure functions at high 
momentum transfer Q for a generic target. First 
we discuss the structure functions in the QCD 
leading order in strong coupling constant (LO) 
in terms of parton distributions (PDFs) focus- 
ing on the analysis of contributions from differ- 



ent isospin states. We point out simple relations 
(see Eas. (|lU|l ) emerging in this approximation. 
Then we discuss perturbative QCD corrections 
to these relations. In Sect.|21we discuss nuclear 
corrections to parton distributions with different 
isospin. The results arc applied to compute the 
non-isoscalarity correction to the ratio /F^ 
for iron. 

2. NC/CC ratios for neutrino scattering 

For an isoscalar target (e.g. the isoscalar com- 
bination of proton and neutron, or for deuterium) 
a relation between neutrino-antineutrino asym- 
metries in the NC and CC deep-inelastic (DIS) 
cross sections was derived long ago by Paschos 
and Wolfcnstein |S] 

R - ^ _ p - - sm Ow, (Ij 
°cc "cc ^ 

where dw is the Weinberg mixing angle. The cor- 
responding relation was also derived for the C- 
even combinations of cross sections (H] (Llewellyn- 
Smith relationship) 

K ■ — ^ - --sm fyvi/ + -^sm ti^ . (2j 
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It should be remarked that in the deriva- 
tion of the Paschos- Wolfenstein (PW) and the 
Llcwcllyn-Smith (LS) relationships the contribu- 
tions from the s and c quarks were neglected. 
Furthermore, the derivation of the LS relation- 
ship holds for high momentum transfer and 
this relationship should be corrected for pertur- 
bative and non-perturbative QCD effects even in 
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an ideal world with only u and d quarks. The PW 
relationship is, however, more general. If only the 
contributions due to light quarks are taken into 
account, the PW relationship is a direct result 
of the isospin symmetry. This ensures that var- 
ious strong interaction effects, including nuclear 
effects, cancel out in R~ for an isoscalar target 
thus making Eq.QJ a good tool for the measure- 
ment of the mixing angle in neutrino scattering. 

If s and c quarks are taken into account the C- 
even ratio i?+ involves contributions due to s -I- s 
and c+c distributions, while the C-odd is cor- 
rected by s — s and c — c asymmetries in the tar- 
get (for a discussion of possible asymmetry in the 
strange sea and the magnitude of this correction 
to R~ see |7I8I9I10) 1 It should be remarked that 
relations Q and (01 are also violated by isospin- 
violating effects in PDFs (for a recent discussion 
of this effect see 

The targets used in neutrino experiments are 
usually heavy nuclei, such as iron in NuTeV ex- 
periment Heavy nuclei typically have an ex- 
cess of neutrons over protons and, therefore, are 
not isoscalar targets. For a non-isoscalar target 
relations ^ and (j2Jl are violated by contributions 
due to isovector component of PDFs. 

In this paper we consider the ratios R^ for 
inclusive differential cross sections for a generic 
target as a function of Bjorken x and the energy 
transfere in units of beam energy E, y = qo/E. 
In terms of the structure functions F2, F^, and 
Fl we have 



R- = 



xFiY^/2 



xF^'Y- 



AF^Y+ 



R+ 



{FiY+ + ^ff y^)/2 
F^Y+ + ^F^y^ + AxF^Y. 



(3) 
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where Y± = ^[1 ± (1 — y)^], the superscript Z 
and W label the NC and CC neutrino struc- 
ture functions, respectively, F^ ~ ^(-F^f + 



and AF^ 



F2 ) and similar defini- 



tions for and the longitudinal structure func- 
tion Fl (in Eqs.© and we neglect the factor 
(1 + QVA/^)V(1 + Q'^jMlf arising due to the 
ratio of propagators of W and Z whose effect is 
small for feasible Q^). 



2.1. Relations between CC and NC struc- 
ture functions 

We now address relations between the NC and 
CC structure functions. We assume that is 
high enough to apply the leading twist (LT) QCD 
approximation. In this approximation the NC 
and CC structure functions are given in terms 
of PDFs. In order to facilitate discussion of 
isospin effects, we consider the isoscalar, q^ix) = 
M(x)-f (i(a;), and the isovector, q\(x) = u(x)—d(x\ 
quark distributions (for simplicity of notations, 
we suppress the explicit notation for the de- 
pendence of parton distributions). We also in- 
troduce quark distributions with definite isospin 
/ = 0, 1 and C parity 



F 



(x) = xqi{x) ± xqi{x), 



(5) 



where q is antiquark distribution. We first con- 
sider the CC structure functions. In terms of the 
functions F'^^^'-^^ the structure functions can be 
written as follows (in QCD leading order) 
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where for C-even and C-odd distributions of 
strange {I = s) and charmed (/ = c) quarks we 
use notations similar to Eq.®. 

For the neutrino NC scattering the LO struc- 
ture functions can be written as 



= C°F^(°'+) + 
C^F^'^+^ + C!2F^'-+\ 
xFi ^C^F^^--^ +C^F^^--^ + 

The coeficients and Cg in Eqs.Q are 



C^ = l-2 sin^ + f sin'' dw, 
= -| sin^ 6*^(1-2 sin^ Ow), 
C| = 1 - I sin^ Ow + I sin^ Ow, 



(7) 
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Using equations © to Q we arrive at the fol- 
lowing relations between CC and NC structure 
functions 



xFi 



(10a) 
(10b) 



These relations have been derived in the LO 
approximation. It is important to study QCD 
corrections to these relations. This is most con- 
veniently done in the DIS scheme in which 
only Fl and F3 structure functions change. In 
this scheme perturbative corrections to Eqs.l^UJ 
are determined by the corresponding correction 
to the structure function ^3. Let us denote 
for each of the structure functions i = 2,3, L 
Fi = Fl^ + SF^, where are the LO structure 
functions by Eqs.® and Q and SFi arc pertur- 
bative series in as which can be written as convo- 
lutions of coefficient functions with LO structure 
functions. For the structure function Ft, 



SF. - 



'-K,iz,asiQ'))Ft''{-,Q% 



(11) 



where K3, is the corresponding coefficient func- 
tion. In the DIS scheme SF2'^ = to all order 
in perturbation theory. Therefore, the general- 
ization of Eq. (|10a|) can be written as 

Fi - C'.fY = -C]AxF3^(^o)(x,Q2), (12a) 

where from the left we have the full (ag-correct- 
ed) structure functions and from the right the 
expression by Eq.® should be used, that is indi- 
cated by the superscript LO. 

The generalization of Ea. ljlOb|) is somewhat 
more complex. We use Eas- Hllfl and Ijl0b|) for 
the LO structure functions in order to find 



xF? - C^xF"^ = -Cl ( AR^(^°' 



[' dzK,{z,as)AF^^'^'=>\-,Q 



(12b) 



where, similar to Ea. H12a|l . from the left we have 
the full structure functions and from the right the 
parton-model expression should be used. 

Equations (|10a|) and (jlObp [as well as as-cor- 
rected Eas. H12a|l and ljl2b|) ] hold for a generic tar- 
get, not necessarily isoscalar. We note that strong 
interaction corrections come through the isosvec- 
tor PDF and/or s or c quark distribution. For an 
isoscalar target the isovector distributions vanish 
F^^'^) = 0, provided that the isospin symmetry 
is exact. If we further neglect the contributions 
from s and c quarks, then the asymmetries Ai^2,3 
vanish and Eas. (|12(l reduce to the LS and PW 
relationships for the structure functions 



Fi/Fi^ = Cl 



Fi/F^- 



(13a) 
(13b) 



Relationships and |(5J) follow if we additionaly 
neglect contributions from the longitudinal struc- 
ture function. It must be remarked, however, 
that the accuracy of the relations in the C-odd 
channel is higher, because they involve asymme- 
tries s — s and c — c, while the C-even relations 
are strongly corrected by strange quark effect at 
small X. Note also that if s = s and c = c (which 
is widely employed assumption in PDF analyses) 
then AF]^ = and Eas. lfT3bjl and ||lj hold to all 
orders in as for the isoscalar target, as it follows 
from Ea. (|12b|l . We also note that 0{as) correc- 
tions for the NC/CC ratios of (anti)neutrino cross 
sections for the isoscalar target were recently dis- 
cussed in Ref . ^21 j in which it was argued that the 
as correction is suppressed by a factor of sin^ 9w- 
For the full analysis one has to take into account 
electro- weak corrections as well |14I15| . 

If the target is not isoscalar then isovector 
quark distributions F^^'^^ may be ffiiite, which 
is the case for complex nuclei. The isovector 
quark distributions violate the LS and PW re- 
lationships for structure functions and cross sec- 
tions. In order to address this effect, in Sect. 01 
we discuss nuclear effects for the isoscalar and 
isovector quark distributions. 

3. Nuclear PDFs with isospin and 1 

Complex nuclei, such as iron, have unequal 
number of neutrons {N) and protons {Z), and the 
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isovector quark distribution is finite in such nu- 
clei. In order to quantitatively understand this ef- 
fect, we denote qa/T ^ the distribution of quarks 
of type a in a target T. If x is large enough to ne- 
glect coherent nuclear shadowing effect, the Icp- 
ton scattering off a nucleus can be well approx- 
imated as incoherent scattering off bound pro- 
tons and neutrons (for a recent review of nuc- 
lear deeply inelastic scattering see CHI)- A widely 
used approximation is to neglect final state inter- 
actions of produced hadronic states with recoiling 
nucleus. In this approximation the nuclear parton 
distribuitions (nPDF) can be written as 



(la/A - {(la/p) + {<la/n)„ 



(14) 



where the two terms in the right side are the 
quark distributions in bound protons and neu- 
trons averaged with the proton and neutron 
nuclear spectral functions, respectively. Similar 
equations can also be written for antiquark distri- 
butions. The explicit expression for the averaging 
in Ea. ((Tl|l is (for derivation and more details see 
|lVl2()l28n 



x{qa/p)p= I ded^kVp{e,k)il + k,/M) 
i/p{ 



2k -q l + {e + k,)/M' 



(15) 
(16) 



The integration in Ea. (|15|l is taken over the en- 
ergy and momentum of bound protons (we sep- 
arate the nucleon mass M from the nucleon en- 
ergy ko = M + e). The quantity Vpie, k) is the 
nuclear spectral function which describes the dis- 
tribution of bound protons over the energy and 
momentum in nuclear ground state. The spectral 
functions Vp and Vn are normalized to the proton 
and neutron number, respectively. In Ea. (|15|l . 
the z-axis is chosen in the direction opposite to 
the momentum transfer q ~ {qo,Oj_, ~\q\), and 
x' is the Bjorken variable of the bound proton 
with four-momentum k. Since bound nucleons 
are off-mass-shell particles their quark distribu- 
tions generally depend on nucleon virtuality k^ 
as additional variable (off-shell effect in structure 
functions is discussed in terms of different ap- 



proaches in |18I19I2()I21I22I230 . Equation sim- 
ilar to Ea. (|15|l also holds for neutrons with the 
obvious replacement of the spectral function and 
quark distributions. 

Equations (|14|1 and H15() account for nuclear 
binding and Fermi motion effects in nPDF (for 
this effect we will use the abbreviation FMB). For 
the isoscalar and isovector nuclear parton distri- 
butions we obtain from Ea. H14|l 



qi/A = {Qi/p)i , 



(17a) 
(17b) 



where the averaging is respectively performed 
with isoscalar and isovector spectral functions, 
Vo = 'Pp+ Vn and Vi=Vp- Vn- 

The isoscalar and isovector spectral functions 
Vo and Vi are very different in complex nuclei. In 
an isoscalar nucleus with equal number of protons 
and neutrons one generally assumes vanishing Vi 
and nuclear effects are dominated by the isoscalar 
spectral function. In a nuclear mean-field model, 
in which a nucleus is viewed as Fermi gas of nucle- 
ons bound to self-consistent mean field, the spec- 
tral function can be calculated as 



VMF{e,p)= ^ nx\(t)\{p)\^ 5{e - £x), 



(18) 



A<A^ 



where 4'\{p) is the wave function of the single- 
particle level A in nuclear mean field and n\ is 
the number of nucleons on this level. The sum in 
Ea. (|18|l runs over occupied single-particle levels 
with energies below the Fermi level \f- Equa- 
tion H18() gives a good approximation to nuc- 
lear spectral function in the vicinity of the Fermi 
level, where the excitation energies of the resid- 
ual nucleus are small. As the separation energy 
e| becomes higher, the mean-field approximation 
becomes less accurate. High-energy and high- 
momentum component of nuclear spectrum can 
not be described by the mean-field model and 
driven by correlation effects in nuclear ground 
state as witnessed by numerous studies (see, 
e.g., |24I25| ). We denote this contribution to the 
spectral function as 7'cor(£,p)- 

For a generic nucleus the spectral function Vi 
determines the isovector nucleon distribution. We 
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now argue that the strength of Vi for complex nu- 
clei is peaked about the Fermi surface. It is rea- 
sonable to assume that Vcor is mainly isoscalar 
and neglect its contribution to Vi. Then Vi is 
determined by the difference of the proton and 
neutron mean- field spectral functions. If we fur- 
ther neglect small differences between the energy 
levels of protons and neutrons then Vi will be 
determined by the difference in the level occupa- 
tion numbers n\ for protons and neutrons. Be- 
cause of Pauli principle, an additional particle can 
join a Fermi system only on an unoccupied level. 
In a complex nucleus all but the Fermi level are 
usually occupied (the Fermi level has a large de- 
generacy factor). Therefore, Vi is determined by 
the contribution from the Fermi level and we can 
write 

Vi^iZ -N)\Mp)\^Sie"eF), (19) 

where ep and (j)F are the energy and the wave 
function of the Fermi level. 

Figure n illustrates isospin dependence of nuc- 
lear effects for valence quark distributions. Shown 
are the ratios = '"V(^^i"'"') and = 
F^J;'^^/[iZ-N)F^^-~^] calculated for the nucleus 
^^Fe26- For the isoscalar nuclear spectral func- 
tion we use the model spectral function of Ref. 
(see also [221 for a different way of fixing param- 
eters of the spectral function) which takes into 
account both the mean-field and correlated con- 
tributions. The isovector spectral function was 
calculated by Ea. H19|l using the Fermi gas model 
in which the wave function of the Fermi level 

I</'f(p)P ~ ^{p ^ Pf)/{^'^Pf), where pp is the 
Fermi momentum. In numerical estimates we use 
£p = -10 MeV and = 260 MeV for the Fermi 
energy and momentum. 

As an application of the present analysis we cal- 
culate the isovector correction to relation (|13b|l . 
Using Ea. l|12b|) and neglecting as correction we 
obtain for the ratio of NC and CC nuclear struc- 
ture functions 

= 1 - sin^ 9w (2 - ISNRf^.ix)) , (20) 

where SN = {N — Z)/A is fractional neutron 
excess in a nucleus and R^/q is the ratio of re- 
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Figure 1. Fermi motion and nuclear binding ef- 
fects for isoscalar and isovector quark distribu- 
tions in ^^Fe26 calculated in impulse approxima- 
tion for Q2 = i5GeV^ 



duced nuclear isovector and isoscalar quark distri- 
butions, [F^'~^/{Z-N)]/{F^'~^/A). The latter 
can also be written in terms of the corresponding 
ratio for the proton R^^q = Fp^'"'' / Fp^'~'^ and the 
ratios i?i and Rq as -Rf/g ~ ^i/o^t/^o- Figure 
2 shows the x dependence of Rf/Q for the iron 
nucleus and the ratio ^^^q for the proton. 

The parameter controlling the isovector correc- 
tion in Ea. H20() is 5N. This correction is addition- 
ally suppressed by the factor of sin^ 9w ■ Given 
that SN w 0.07 for iron we observe that the 
isovector correction in Ea. (|20|l is not large. In 
particular, in the valence quark region of 2; ~ 0.3 
the magnitude of the correction in Ea. (|20|l is 
about 1% of the PW value 1-2 sin^ Ow This cor- 
rection rises with x, as is clear from Fig.[5| How- 
ever, its magnitude is only about 2% at x = 0.7. 

We note that the non-isoscalarity correction to 
R~ for the total cross sections was recently dis- 
cussed in |27I28| . It was observed that the proper 
treatment of the FMB effect leads to about 6% 
enhancement of the neutron excess correction for 
iron nucleus. We also remark that a subtle cancel- 
lation of perturbative QCD effects in the isovector 
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Figure 2. The ratio of isovector and isoscalar 
quark distributions calculated for the proton and 
iron at = 15 GcV^ 



correction for R~ was found in |28j . 

In summary, wc discussed relations between 
NC and CC structure functions for a generic 
nuclear target. We argued that the relation 
^3/^3^ ~ ^ ^ 2sm^9w survives strong inter- 
action corrections if the strange sea in the target 
is symmetric and the target is isoscalar. The cor- 
responding relation for F2 is less accurate and 
affected by strange quark effect at small x. We 
also discussed how the parton distributions with 
isospin and 1 change in complex nuclei and ap- 
plied the results to compute the non-isoscalarity 
correction to the ratio /F^ . 
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